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ABSTRACT 

With a non local shell model of magnetohydrodynamic turbulence we investi- 
gate numerically the turbulent dynamo action for low and high magnetic Prandtl 
numbers (Pm). The results obtained in the kinematic regime and along the way 
to dynamo saturation are understood in terms of a phenomenological approach 
based on the local {Pm ^ 1) or non local {Pm ^ 1) nature of the energy trans- 
fers. In both cases the magnetic energy grows at small scale and saturates as an 
inverse " cascade ". 

Subject headings: MHD, turbulence, methods: numerical, magnetic fields, plas- 
mas 



1. Introduction 

Dynamo action is generally believed to be at the origin of the magnetic field in most 
astrophysical objects. The conducting fluid which produces the magnetic field, is generally 
strongly turbulent {Re ^ 1) with an extended inertial range. In addition as it is electrically 
conducting, the fluid is characterized by the magnetic Prandtl number Pm = v/rj, where v 
and rj are the fluid viscosity and magnetic diffusivity. In a liquid metal as in planetary cores 
or stellar convective zones Pm -C 1, while in the warm interstellar medium, coronal and 



cluster plasmas Pm ^ 1 fISchekochihin etal. 20021 ). In direct numerical simulations both 
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conditions, Re ^ 1 and Pm much different from unity, are still out of reach of the present 
day computers. An alternative is to use shell models. 

Hydrodynamical shell models are a rough approximation of the Navier-Stokes equations 
formulated on a discrete set of real wave numbers kn = A", corresponding to a Fourier 
space divided into shells of logarithmic width A. They involve (complex) scalar quantities 
reminiscent of the velocity Fourier components of an isotropic flow. The model satisfies the 
conservation of both kinetic energy and kinetic helicity when the viscosity is set to zero. It 
leads to the resolution of a system of ordinary differential equations, requiring much less 
computing power than the direct numerical simulations of the original equations. It is then 
possible to run a shell model for realistic values of the viscosity such as u = 10~^m^ ■ (for 
a review on shell models, see Biferale (2003) and references therein). 

The results obtained with shell mo dels are in general ag reement with the Kolmogorov 



phenomenology including intermittency (ILeveque fc She 19971 ) . Firstly at any scale k ^ lying 



in the inertial range, the flux rate of kinetic energy is equal to the injection rate of kinetic 
energy at the forcing scale. This writes ku^{k) ~ e, leading to u(/c) ~ e^^^k^^/^ and to the 
Kolmogorov spectrum E{k) ~ e^^^k^^^^. Secondly, the inertial range is found to extend to the 
viscous scale k~^ which can be evaluated saying that at this scale the turn over time k~^u^^ 
compares with the viscous time z/~^fc~^, leading to k^j Thirdly at low viscosity 

some deviation of the Kolmogorov power scaling is found due to intermittency. The power 
scaling of the statistical moments are given by < \u{k)\ >p~ with scaling exponents 
deviating from Kolmogrov's mean field theory C,p = p/3. The value of these anomalous 
scaling exponents giy en by the shell models compare well with those found in experiments 



( ILeveque fc She 19971 ). Therefore shell models appear to be a useful tool to study fully 



developed turbulence at hi gh Reynolds nurnbers. They are also used in their MHD version t o 



tackle astrophysical issues flFrick etal. 200fil : ISuchlin fc VeUi 20071 : lOaltier fc Buchlin 20071 ). 



In a previous paper (jPlunian fc Stepanov 20071 ). denoted after by PS07, we introduced 



a non local shell model which is not only in agreement with the previous features of turbu- 
lence but in addition gives an appropriate slope of the infrared spectrum in freely decaying 
turbulence (which is not the case of the local shell models described in Biferale 2003). The 
magnetohydrodynamic (MHD) version of this model, also introduced in PS07, permits to 
study turbulent dynamo action for arbitrary low or high values of Pm. In PS07 we calcu- 
lated the energy transfer functions of the MHD system in a (statistically stationary) satu- 
rated state. We found that for Pm < 1 the energy transfers are mainly local, eventually 
strengthening our previous results obtained with a local shell model of MHD turbulence 



( jStepanov &: Plunian 20061 ). For Pm ^ 1 the dominant transfers are also mainly local ex- 



cept the ones from the flow scales lying in the inertial range to the magnetic scales smaller 
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than the viscous scale. In that case the use of a non local model is definitely necessary. 

Here our goal is to describe the transient properties of the turbulent dynamo, namely 
the kinematic regime during which the magnetic energy grows exponentially, and the route 
to saturation, starting when the kinematic regime stops and ending when a (statistically 
stationary) saturated state is reached. For that it is not sufficient to make time averages as 
we did in PS07 when studying the saturated dynamo states. Indeed the transient properties 
we are interested in are, by definition, not statistically stationary. Then instead of time 
averages we have made ensemble averages over M = 12 ■ 10^ independent realisations. The 
initial conditions have the same energy but random phases. In addition the forcing has also a 
random phase in time (therefore different from one realisation to the other). Then the results 
between two realisations are different. Every to = 10~^ unit of time we freeze the calculation. 
For each shell we average the magnetic and kinetic energy over the M realisations without 
modifying the phases of the magnetic, kinetic and forcing quantities in order to keep the 
randomness between every realisations. All the calculation were done for = 60 shells, 
corresponding to more than 12 decades. Before presenting the results (sections |3] and H] 
for low and high Pm) we briefly describe in the next section the shell model (previously 
introduced in PS07). 



2. The non local shell model 

The model is deflned by the following set of equations 

Un = tK[QniU,U,a)-QniB,B,a)]-iyklUn + F^, (1) 
Bn = ikn[QniU,B,b)-QniB,U,b)]-r]klBn, (2) 

where 

N 

Qni^i ^) c) = Tjn[Cj^X^^^n+m+l + C^Xn-mXri+l + (^^n-m-l^n-l\ (3) 

m=l 

represents the non linear transfer rates and F„ the turbulence forcing applied at shell n = 0. 
As explained in PS07 an optimum shell spacing is the golden number A = (1 + \fE)/2. 



For A^ = 1 in ([3]), we recognize the local Sabra model (ILVov etal. 19981 ). The additional 
non-local interactions for N > 2 correspond to all other possible triad interactions except 
the ones involving two identical scales. Expressions for the kinetic energy and helicity Eu 
and Hu, magnetic energy and helicity Eb and Hb, and cross helicity He are given by 

Eu = J2Mn), Eu{n) = \\Un\\ Hu = J^^^in), Hu{n) = 1(-1)"A;„|[/„|2 (4) 
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Eb = J2EB{n), Esin) = ^\B^\', Hb = Y,HB{n), HB{n) = l(-l)"A;-i|5„|2 (5) 

n n 

He = J^Hcin), Hc{n) = \{UnBl + B.^U:). (6) 

n 

In the inviscid and non- resistive limit (z/ = r/ = 0), the total energy E = Ejj + Eb-, 

magnetic helicity and cross helicity must be conserved {E = Hb = He = 0). This implies 
the following expression for the coefficients and 6^: 

„1 _ L , L „2 _ -Wi-(-l)'" „3 _ fcm - (-l)"" 

ITT- ' ""m+l k m 



h]n = {-ir^' hl = l hl = -l. (7) 



m 



In the case of pure hydrodynamic turbulence (without magnetic field), the coefficients a 
derived from the kinetic energy and helicity conservations {Eu = Hjj = 0), would lead to 
the same expression as ([7]). The coefficients are free parameters depending on m only, 
that we choose of the form = k^_i/X{X + 1). Further characteristics of the model and 
results can be found in PS07. In particular the role of the non locality parameter a has 
been investigated. It is indeed the only free parameter left in the model that can not be 
theoretically constrained. 

In PS07 we found that taking a = —5/2 permits to describe accurately the infrared 
hydrodynamic spectrum in freely decaying turbulence. Keeping this value of a in the MHD 
system for Pm < 1, we found that the dominant energy transfers are mainly local. 

On the other hand, for Pm ^ 1, in order to describe properly the non local energy 
transfer from the flow scales lying in the inertial range to the magnetic scales smaller than 
the viscous scale, the results obtained in PS07 suggest to take a parameter a corresponding 
to stronger non local transfers than those obtained for a = —5/2. A phenomenological 
justiflcation of the choice a = —1 will be given below. It can be also understood as a way to 
mimic in our (isotropic) model the strongly anisotropic small scale magnetic turbulence. It 
is worth mentioning here that taking a = — 1 does not affect too much the energy transfers 
involving the magnetic scales lying in the inertial range, which remain mainly local. 



3. Low Pm dynamo action 

3.1. Kinematic regime 

At low Pm, as shown in PS07, the energy transfers responsible for the generation of the 
magnetic from the kinetic energy are mainly local and the flow scales producing the largest 



- 5 - 



magnetic growthrate lie in the inertial range. A phenomenological way to understand it is 
to write the induction equation at wave number k, in the form 

b{k) = k U{k)b{k) - 7]k%{k) (8) 

leading to the following growthrate 

-f{k) = ku{k) -r]k'^. (9) 

Assuming a Kolmogorov spectrum in the flow inertial range, u{k) = e^^^k"^^^, we find 
7(fc) = e^/^k"^/^ — rfk'^ which is plotted in figure [H In the kinematic regime due to the linearity 
of the induction equation all magnetic scales grow with the same growthrate 7fcj„ = max7(fc). 
It occurs at k^in = 3~^^'^€^^'^r]~^^'^ with the value jkin = 2 / {3\/3)e^^'^r]~^/'^ . The resistive scale 
k~^ is the scale for which the shear ku{k) is compensated by the resistive dissipation r^/c^, or 
equivalently it is the scale for which 7(fc) = 0. We have k^ = €^^^r]~^^^ and ki^in = 0.438 fc^. 

Defining a magnetic Reynolds number at wave number k by Rm{k) = k^^u{k)/r], 
we have Rm{k) = e^^^k~'^^^ /r] in the inertial range and consequently Rm{kn) = 1 and 
Rm^kkin) = 3. Comparatively the magnetic Reynolds number at the forcing scale kp^ 
is much larger, scaling as Rm{kp) ~ Rm{kkin)/r]. Surprisingly the kinematic growthrate 
is given by the small scale kj^^^ for which Rm{k) is of order unity and not by the forc- 
ing scale kp^ for which Rm{k) is much larger. In addition kkin and jkin are viscosity 
independent in the limit of low Pm, suggesting that additional scales in the kinetic in- 
ertial range do not change the result. This is also consistent with the fact that the dy- 
namo threshold should not depend on Pm in the limit of low Pm as shown previously 



flPontv etal. 20051 : IStepanov fc Plunian 2009 : ISchekochihin etal. 20071 1. 



Coming back to the shell model, as the induction is linear in B = {bi,b2, ■■■ ,bi\i) where 
is the maximum number of shells, we end up with a system of the form dB/dt = M. ■ B 
where is a A^ by A^ matrix depending on the flow U = {ui,U2, ■ ■ ■ ,un)- During the 
kinematic regime the flow U is statistically stationary (the kinematic regime being defined 
as long as the Lorentz forces are negligible), implying that the system dB/dt = M. ■ B can 
be solved as an eigenvalue problem. Among the A^ eigenvalues of the matrix A4, the largest 
one should correspond to 7fcj„ (provided that our phenomenological approach is correct). In 
addition the corresponding eigenvector should be peaked at kkin- 

In figure[2^ l""!^)!^ (black curves) and |6(fc)p (red curves) are plotted versus k at different 
times, for u = 10~^ and r] = 10~^. The kinetic energy satisfies the Kolmogorov scaling 



u 



2 



(some deviation occurs due to intermittency as shown in PS07). In the 
kinematic regime the magnetic spectrum (red curves in [2^) is peaked at a scale in agreement 
with kkin and grows exponentially with the same growthrate at all scales as shown in figure 
[2}1, with a growthrate in agreement with •jkin- 
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Fig. 1. — Growthrate spectrum ^{k) produced by u{k) lying in the inertial range. 
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Fig. 2. — Low magnetic Prandtl number regime Pm -C 1. In all figures except {h), u = 10"'' 
and 1] = 10^"^. In (a) and (c), \u{k,t)\'^ (black curves) and |6(/c,t)p (red curves) are plotted 
versus k at different times. In (d) \b{k,t)\^ is plotted versus time for several values of k. 
The thickest (and darkest) curve corresponds to the largest scale. The magnetic energy 
growthrate "^{k^t) is plotted versus k in (e) and versus t in (f). In (e) the curves evolve from 
top to bottom along time. In (f) the color scale is identical to (d). In (b) the growthrate 
during the kinematic regime is plotted versus rj for several viscosities. 
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The kinematic growthrate has been calculated for other values of Pm < 1. It is plotted 
versus t] in figure [2)d for three viscosities (z/ = 10~^, 10~^, 10^'^). It is clear from the over- 
lapping curves that the growthrate does not depend on the viscosity. In addition the curves 
show a clear power dependence in t]. The above phenomenology tells us that 'jkin ~ with 
r = —0.5. Instead we find from figure [2)d that r = —0.46. This apparent discrepancy can be 
explained from the fact that, highly intermittent, we have m„ oc /c„^ with ( = 0.369 

(instead of the Kolmogorov scaling ( = 1/3). This new scaling implies ~ 7^-1/(1+0 ^^d 
7fcm ~ r7-(i-^)/(i+«), leading to r ^ -0.46. 



3.2. Dynamic regime 

The kinematic regime ends when the saturation starts. Here this corresponds to the 
time when the magnetic energy at wave number kkm becomes comparable to the kinetic one 
b{kkin) ~ u{kkin)- Then due to the Lorentz forces the magnetic energy at that scale saturates 
and its growthrate decreases to zero. Meanwhile, the magnetic energy of the next smaller 
wave number (fcfcin/ A in the shell model) becomes the one with the largest growthrate. Then 
the wave numbers k < kkin/^ continue to grow with the new growthrate 'y{kkin/^) until 
h{kkin/X) saturates b{kkin/X) ~ u{kkin/X). Then the next smaller wave number kkin/^^ 
becomes the one with the largest growthrate and successively all the next smaller wave 
numbers change their growthrate in accordance with the growthrate spectrum 7(/c) given in 
figure [H This saturation mechanism lasts until the smallest wave number is saturated. Then 
equipartition between both kinetic and magnetic energies is obtained at all scales. 

Such a scenario is quantitatively illustrated by the mode crossing in figure [2](d) where the 
magnetic energy is plotted versus time for several wave numbers. The way the magnetic wave 
numbers get their energy is reversed when compared to the direct hydrodynamic cascade, 
be cause it starts at kk in towards smaller wave numbers. We speak of an inverse " cascade 



" (jPouquet etal. 19761 ). It is also illustrated in figure [2](e) where the magnetic growthrate 
t) /b{k, t) is plotted versus k at different times. In the kinematic regime the growthrate is 
the same at all wave numbers (top black curve). Then as time goes, the plateau of constant 
growthrate decreases in both intensity and range of involved wave numbers. 

An interesting feature of the route to saturation is shown in figure [2](c). The injection 
rate of kinetic energy e partly dissipates at the resistive scale k~^. A direct consequence is 
that the viscous scale k~^ shifts to larger values during the saturation regime. 

Finally in figure [2](f) the growthrate b{k, t)/b{k, t) is plotted versus time for each scale k. 
A signature of the saturation scenario described above is the time evolution of the growthrate 
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at the largest scale (black curve). We shall compare it to the following phenomenological 
description which is sketched schematically in figure [31 

At time t the magnetic field at wave number k reaches its saturation value b{k, t) = u{k). 
The left part of the magnetic spectrum scaling as A;", we have b{k + dk, t) = b{k, t){l + dk/k)"'. 
Then for t <t' < t+dt the magnetic field at scale k + dk grows exponentially as b{k+dk, t') = 
b{k + dk,t) exp{'j{k + dk) ■ {t' — t)). At time t' = t + dt it reaches its saturation value 
b{k + dk,t + dt) = u{k + dk). For a Kolmogorov turbulence u{k + dk) = u{k){l + dk/k)~^/^ 
and assuming that dk"^ + dt^ <^ 1, we finally find 

From ([9]) and for k < kkin-, a rough approximation of 7(/c) is 7(A;) ~ 'Jkinik/kkm)'^^^ leading 
to 

dk 'Jkin 



dt {a + l/3)ky.^ 



(11) 

Defining to as the time at which the saturation starts (at k = kkm) , we find that the magnetic 
scale which saturates at subsequent time t > to satisfies 



-3/2 



kkin \ (1 + 3a) / 
The growthrate of the largest magnetic scale then satisfies 

This is the signature of the saturation scenario that we want to compare to the shell model 
results. In the shell model we have a ~ 2 (corresponding to an " infrared" spectrum b'^{k) ~ 
/c^ as found e.g. in Christensson etal. 2001). Therefore we expect the inverse of the magnetic 
energy growthrate in shell n = to satisfy (27)"-^ = (27fcj„)~^ + (t ~ ^o)/7. In addition 
intermittency again changes slightly the result. Indeed for u ~ k~'^, we find that 

7-'(t)=7,-i + ^^^(t-to) (14) 

leading to a slope d{2'y)^^/dt = 0.133 instead of 1/7 (without intermittency). In figure HI 
(27)-^ is plotted versus time. We find indeed a good agreement with the previous expecta- 
tion, at least for 9 < t < 11. At subsequent times the curve diverges because all magnetic 
wave numbers are saturated expect the 2 smallest ones {n = 0, 1). As the turbulence forcing 
is applied at these wave numbers, they do not belong to the inertial range and then the 
previous picture fails. Therefore the saturation scenario works until the saturation scale 
becomes close to the forcing scale. 
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4. High Pm dynamo action 

4.1. Kinematic regime 

At high Pm, the resistive scale is much smaller than the viscous scale k~^. The 
flow scale which is responsible for the magnetic growth in the kinematic regime is the one 
which has the highest shear ku{k). Assuming a Kolmogorov spectrum in the inertial range, 
we have ku{k) = e^^^k"^^^ which is maximum at the viscous scale k~^ (at smaller scales the 
viscous dissipative range starts and u{ k) falls down). Then the magnetic growthrate in the 



kinematic regime can be estimated as (jSchekochihin etal. 20021 : ISchekochihin etal. 20041 ) 



7 = k,u, = {e/uy/\ (15) 

We note that 7 depends on the viscosity but not on the diffusivity. The corresponding 
phenomenologic induction equation is 

b{k) = k^Uyh{k) - r]k%{k) for k > k^ (16) 



The resistive scale k~^ is the scale for which the shear k^^Uy compensates for the resistive 
dissipation 77 A;^. As at the viscous scale the shear also compensates for the viscous dissipation, 
we have rfk'^ = uk^, leading to k^^/k^ = Pm}^"^. 

From f|T6|) we expect a flat spectrum of |&(/c)P for k^ < k < krj. In addition the "infrared" 
spectrum k < k^ is enslaved to the growth of the largest wave numbers k > k^. 

Comparing ([2]) with (1161) we can show that the only way to account for the non local 
term ki,Uyh{k) is to take a = — 1 in the shell model. In flgure[5](a), |'u(A;)P (black) and |&(/i;)P 
(red) are plotted versus k at different times, for v = 10^^ and 77 = 10"*^. In the kinematic 
regime, the magnetic spectrum is not completely flat for k^ < k < k^. This is probably an 
effect of the additional non local energy transfers involving kinetic scales larger than k~^. 

The kinematic growthrate that we flnd agrees with f|T5l) . It has been calculated for 
other values of Pm > 1. It is plotted versus v in flgure[5t^b) for four values of the magnetic 
Prandtl number {Pm = 10,10^,10^,10^). For increasing values of Pm the curves shift 
from bottom to top and reach an asymptotic limit showing that the growthrate becomes 
?7-independent at high Pm. In addition the curves show a clear power scaling in u. The 
above phenomenology tells us that '-fkm ~ with r = —0.5. Instead we flnd from flgure[5](b) 
that r = —0.46. As in the low Pm case, this apparent discrepancy results from the fact that 
the velocity Un is intermittent. Then for oc k~'^ with ( = 0.369, we flnd k^ ~ jy-^/i^+O 
and 'jkin ~ i^-^'^-O/i^+O ^ leading to r ^ —0.46. 
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Fig. 5. — High magnetic Prandtl number regime Pm ^ 1. In all figures except (h), u = 10~^ 
and rj = 10~^. The caption is the same as in figure [2] except (b) for which the growthrate 
during the kinematic regime is plotted versus u for several diffusivities. 
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4.2. Dynamic regime 

As explained before, during the kinematic regime, the magnetic field b{k) at wave num- 
ber k > receives energy from Uy. It also releases energy to u{k) locally such that the 
Lorentz force partially compensates for the viscous dissipation. This is visible in figure [5]^a) 
where in the viscous range for scales k^ < k < fc^, the slope of becomes less steep as 

soon as the magnetic energy is sufficiently large. 

At wave number k = k^ when \b{ku)\'^ ~ \u{ky)\'^ the saturation starts and again the 
saturation occurs from the smaller to larger scales as shown in figure [5] (c) to (e). Again the 
inverse "cascade" saturation scenario seems to be the relevant one. 

In figure [6], (27)"^ is plotted versus time. For 3.5 <t< 7.7 we find a good agreement 
with the phenomenological prediction derived in section 13.21 As explained previously, at 
subsequent times the curve diverges as the remaining scales to saturate are larger than the 
inertial range. 



Discussion 



In the kinematic regime, we found that the magnetic growthra te satisfies 7 ^ (e/vY^'^ for 
Pm <^ 1 and 7 ~ (e/i/)^/^ for Pm ^ 1 and is therefore always fast (jchildress fc Gilbert 1995 ) 
In both cases this corresponds to a small-scale dynamo (or fiuctuation dynamo in the termi- 
nology of Schekochihin etal. 2007) in opposition to a large-scale (or mean- field dynamo). At 
low Pm the magnetic spectrum is peaked at A; = kkin ~ OAkrj and the energy transfer from 
kinetic to magnetic is mainly local (as shown in PS07 the non local transfers from small 
kinetic to large magnetic scales is less than 20 % than the local transfers). At high Pm 
the magnetic spectrum is almost fiat between the viscous and resistive scale and the energy 
transfer is mainly non local from u{kij) to b{k) with k^^ < k < fc^. In the dynamic regime 
we found an inverse "cascade" mechanism which explains the route to saturation. The rel- 
evance of this mechanism does not depend on Pm. When all magnetic scales are saturated 
the resulting equilibrium between kinetic and magn etic energy is not equipartition with a 
slight excess of magnetic energy (IHaugen etal. 20031 ). We find that the ratio \b{k)/u{k)\ is 
constan t for kp < k < min|fc„, fc^j corresponding to a residual energy with a Kolmogorov 
scaling (iMiiller &: Grappin 20051 ). 



Dynamo action is usually found in natural objects with strong rotation. A transition 
from a Kolmogorov spectrum E(k) ~ k~^^^ for weak rotation to E(k) ~ for strong rota- 



tion is predicted by weak turbulence theory, and observed in experiment s (IBaroud etal. 2002 



Baroud etal. 20031 ) and numerical simulations (jSmith fc Waleffe 19991 ) including shell mod- 
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els ( iHattori etal. 20041 ) . It is then of interest to see how the previous phenomenology is 
changed. At low Pm for u{k) ~ k~'^ with 1/3 < C ^ 1/2? we derive the quantities given in 
table [T] and give the numerical values expected for the three scalings, Kolmogorov, intermit- 
tent and strongly rotating. The change of the power scalings from weak to strong rotation 
should be sufficient to be identified in measurements or numerical calculations. On the other 
hand kkin/k^q and Rm{kkin) do not change much. 



Most of this work was done during the Summer Program on MHD Turbulence at the 
Universite Libre de Bruxelles during July 2007. The organizers D. Carati and B. Knaepen 
are warmly thanked as well as the sponsors of this program. R.Stepanov thanks for financial 
support from RFBR grants (07-01-96007 ural and 07-02-00127) and the Russian Federation 
President grant MK-4338. 2007.1. The simulations were performed on the computer cluster 
of IMM (Ekaterinburg, Russia). 





Phenomenology 


Kolmogorov 


Intermittency 


Strong rotation 


c = 




1/3 


0.369 


0.5 


kkin/krj 


(¥)* 


0.439 


0.431 


0.397 


Rm{kkin) = 


2 


3 


3.17 


4 


kkinj k^ ~ 


^-0.75 


^-0.73 


^-0.66 


k,y ~ 


1 




^-0.73 




Ikin ~ 


1-C 

r] i+c 


^-0.5 


^-0.46 


^-0.33 


'^(27)-^ _ 
dt 


i-C 

a+C 


0.143 


0.133 


0.1 



Table 1: Summary of the numerical values expected from phenomenological predictions at 
low Pm for different velocity spectra u{k) ~ A;"''. 
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